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Abstract 

We present the results of a quenched lattice calculation of the operator ma- 
trix elements relevant for predicting the B s width difference. Our main result is 
(ATb s /Tb s ) = (4.7 ± 1.5 ± 1.6) x 10 -2 , obtained from the ratio of matrix elements 
K(m h ) = (B s \Qs\B s )/(B° s \Q L \B°) = -0.93(3)^. K(m b ) was evaluated from the 

two relevant ^-parameters Sf^(m 6 ) = 0.86(2) ±°;° 2 and B^f(m b ) = 0.91(3)l£$j, 
which we computed in our simulation. 



PACS: 13.75Lb, 11.15.Ha, 12.38.Gc. 



1 Introduction 



In the Standard Model, the width difference (AT Bs /T Bs ) of B s mesons is expected to be 
rather large and within reach for being measured in the near future. Recent experimental 
studies [[[], |2| already provide an interesting bound on this quantity. In particular, in ref. || 
the limit (AT B jT Bs ) < 0.31 (95%C.L.) is quoted Q. 

Theoretically, the prediction of (AT/T) B relies on the use of the operator product ex- 
pansion (OPE), where the large scale is provided by the heavy quark mass ||. All recent 
developments, including the calculation of the next-to-leading order (NLO) perturbative 
QCD corrections, have been discussed in great detail in refs. The theoretical esti- 

mates are in the range (AT Bs /T Bs ) = (5 -f- 15)% and crucially depend on the size of relevant 
hadronic matrix elements which must be computed non-perturbatively. 

In this paper we present a new lattice calculation of the main contribution to (AT Bg /T Bg ). 
On the basis of our results, and using the expressions given below, we predict 

AY 

— - = (4.7 ± 1.5 ± 1.6) x 1(T 2 , (1) 

where the last error is obtained by assuming an uncertainty of 30% on the l/m b corrections. 

We now present the relevant formulae which have been used to get the prediction in 
eq. ([]]). Up and including l/m b corrections, the theoretical expression for (AT B jT Ba ) 
reads [HI 



Ar« G 2 m 2 



B s _ ^F« l b it/ T/ |2 



V cb V cs \\ Bs ^G(z)(Q L (m b )}-G s (z)(Qs(m b )}+5 1/m VT^) , (2) 



where z = m 2 /m 2 and G(z) and Gs(z) are functions which have been computed in per- 
turbation theory at the next-to- leading order (NLO) ||. (Ql,s) — (Bs\Ql,s\^s) are the 
hadronic matrix elements of the renormalized operators relevant at the lowest order in the 
heavy quark expansion 

ql = &y(i-7 5 y ^v(i-7 5 y 

Qs = 6X1 - 75)^^(1-75)^', (3) 



where i and j are colour indices. The last term in eq. (§), 5\/ m y/l — 4z } represents the 
contribution of l/m b corrections. 

In order to reduce the uncertainties of the theoretical predictions, it is convenient to 
consider the ratio of the width and mass differences of the B° s -B° s system 



AY Bs _ 4tt m 2 b 



Am Bs 3 m 2 



VcbV cs 



V ts V tb 



G(-)-G s (z)^ + ~5 1/m ) . (1) 



rj B (m b )S Q (x t ) V (Ql) 



where r) B (m b ) has been computed in perturbation theory to NLO [0 and S (x t ) is the usual 
Inami-Lim function |@|. Note that, to make contact with ref. [HI, in the above formulae we 



1 For this estimate, the average B s decay width was assumed to be the same as for Bd mesons. 
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used the MS coefficient ^(m?,) instead of the standard r\ B of ref. J7|. Consequently, the 
operators Qs and Ql are renormalized in the MS (NDR) scheme. We see from eq. (f|), 
that (Ar^/AmgJ only depends on the ratio of matrix elements, 



TZ(m b ) 



(Q s ) (B° s \Q s (m b )\B° s } 



(Ql) 



{m b )\B°) 



(5) 



which may, in principle, be directly determined on the lattice. The use of TZ is particularly 
convenient because dimensionless quantities are not affected by the uncertainty due to the 
calibration of the lattice spacing. One may also argue that many systematic errors, induced 
by discretization and quenching, cancel in the ratio of two similar amplitudes. 

Finally, eq. (|) allows to express (ATb s /Tb s ) in terms of i) perturbative quantities, 
encoded in the overall factor and in the functions G(z) and Gs{z), ii) a lattice measured 
quantity, TZ and iii) ArriB s which will be hopefully precisely measured in the near future: 



TBArnB^ Kx (G(z)-Gs(z)n(m b )+5 1/r 



where 



K 



4tt ml 
3 m 2 w 



V cb V cs 



v ts v tb 



r) B (m b )S (x t ) 



Waiting the measurement of ArriB s , for which only a lower bound presently exists 
can use a modified version of eq. (](]), namely 



(6) 



(7) 



one 



Ar, 



B s 



T B Am B , 



m Bs 
' ™b, 



where 



| (exp.) 


v ts 




v td 




_ h 








f B 



K ■ (G(z) - G s {z)n{m b ) +6 1/m ) £ 



B, 



B B d 



(9) 



In this way, besides the quantities discussed above, we only use the experimental £> d -meson 
mass difference, which is known with a tiny error 



(Am Ba 



>(exp.) 



0.484(15) ps^ 1 



(10) 



and another ratio of hadronic matrix elements, namely £, which is rather accurately deter- 



mined in lattice simulations 10, 11 



For the following discussion, it is useful to write eq. (|8|) as (note the 7Z(m b ) is negative) 



B 3 



(0.5 ± 0.1) - (13.8 ± 2.8)K(m b ) + (15.7 ± 2.8)^ 



/in 



x 10" 



(11) 



where the three contributions correspond to G(z), —G s (z)1Z{m b ) and 5i/ m in (||), respec- 
tively. 
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The advantage of using eqs. (Q), (||) and (|8|) consists also in the fact that, in order to 
predict (AT Bs /Tb s ), we do not need f Bs , which enters eq. (0) when we express the matrix 
elements in terms of 5-parameters. There is still a considerable uncertainty, indeed, on f Bs , 
which has been evaluated both in the quenched (f Ba = 195 ± 20 MeV) and unquenched 
(f Bs = 245 ± 30 MeV) case, with a sizeable shift between the two central values ]10| . 



Since the "unquenched" results are still in their infancy, however, we think that the large 
quenching effect should not be taken too seriously yet. 

In the numerical evaluation of the width difference, we have used the values and errors 
of the parameters given in Tab. [I]. For the perturbative quantities G(z) and Gs{z), the 
main uncertainty comes from dependence on the renormalization scale, which was varied 
between nib/ '2 and 2nib in ref. || (with nib = 4.8 GeV). For this reason, we found it useless 
to recompute these coefficients with nib = 4.6 GeV (which is very close to the mass used 
in [BJ). Instead, we took as central value the average G s (z) = (G ,I s 1=mb/2 (z)+G' x s 1=2 m "(z))/2 
(see Tab. 1 of that paper), and as error {G^~ mb ^ 2 {z) — G^~ 2 mb (z))/2. The same was made 
in the case of G(z). 

For the hadronic quantities, we have used the following values 

TZ(m b ) = -0.93(3)1°;°? , £ = 1.16(7) , (12) 
lZ(nib) is the main results of this lattice study, while the result for £ has already been 



reported in our previous paper | IT| . Following ref. Q, for the l/m& corrections we get 

5 1/m = -2.0GeV 4 , (13) 



using f Ba = 204 MeV from ref. JTTJ] . This value of 8u m corresponds to 8u m = —0.55 
(obtained with B B (mb) = 0.91) Q. Since in the estimate of 5i/ m , the operator matrix 
elements were computed in the vacuum saturation approximation (VSA), and the radiative 
corrections were not included, we allow it to vary by 30%, i.e. 5u m — —0.55 ± 0.17. In the 
numerical evaluation of the factor K, we have used the pole mass m& = 4.6 GeV derived 
from the MS mass in Tab. |T] at the NLO. Using these numbers, from eq. (^j) we obtain the 
result in eq. (|l|), where the last error comes from the uncertainty on 8\/ m . 

We have an important remark to make. Eq. ([H]) shows explicitly the cancellation 
occurring between the main contribution, proportional to TZ(nib), and the 1/nib correc- 
tions, proportional to Si/ m . Without the latter, we would have found a much larger value, 
(AT Bs /r Bs ) ~ 14 x 10~ 2 . (AT Bs /r Bs ) would remain in the 10% range even with a size- 
able, but smaller, value for the 1/nib term, 8y m = —0.30 for example. This demonstrates 
that, in spite of the progress made in the evaluation of the relevant matrix elements ((Qs) 
and (Ql)), a precise determination of the width difference requires a good control of the 
subleading terms in the l/my, expansion, which is missing to date. 

One can also use eq. (|), and combine it with 

(Am B J (cxp - } > 14.6 ps' 1 , (14) 



2 In ref. 0, we gave B^f(m b ) = 0.92(3)+°;°6- The tiny difference is due to the fact that there, in the 
perturbative evolution, we used a number of flavours n/ = 4 instead of n/ = as in the present paper. 
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Parameter 


Value and error 


m w 


80.41 GeV 


m Bd 


5.279 GeV 


m Bs 


5.369 GeV 




1.460 ± 0.056 ps 


\Vcb\ 


0.0395 ±0.0017 


\Vu\ 


0.0386 ±0.0013 


\Vcs\ 


0.9759 ±0.0005 


\Vtd\ 


0.0080 ± 0.0005 


u s {M z ) 


0.118 ± 0.003 


Am Bd 


0.484 ±0. 015 ps" 1 


m t 


165 ± 5 GeV 


m b 


4.23 ± 0.07 GeV 


m c /m h 


0.29 ± 0.02 




110 ± 20 MeV 


r)B\m h ) 


0.85 ± 0.02 


G(z) 


0.030 ±0.007 


Gs(z) 


0.88 ±0.14 



Table 1: Average and errors of the main parameters. When the error is negligible it has 
been omitted. The heavy-quark masses (m t , mi, andm c ) are the MS masses renormalized at 
their own values, e.g. m t = mf ts (m^ IS ). m s = m^ s (/i = 2 GeV) is the strange quark mass 
renormalized in MS at the scale fi = 2 GeV. Its central value and error are not important 
in the numerical evaluation of ATb s - 



to obtain a lower bound on (ATbJTbJ- Given the large uncertainties, this bound is at 
present rather weak. At the 1-cr level we get 



Ar, 



r, 



> 2.5 x 10~ 2 (68%C.L. 



(15) 



Following ref. 0, from eq. (^) and the limit (ATbJTb s ) < 0.31, we could also obtain an 
upper limit on Am^. In our case this is not very interesting, however, since we find a very 
large upper bound of O(200) ps^ 1 . 

Our prediction in eq. ([[]) is in good agreement with ref. 0. It is instead about a 
factor of two smaller than the result of ref. ||17|| . A detailed comparison of the two lattice 
calculations can be found in sec. ^. 

We stress that the theoretical formulae should be evaluated with hadronic parameters 
computed in a coherent way, within the same lattice calculation, and not from different cal- 
culations (the "Arlequin" procedure according to ref. [|12j ) , since their values and errors are 



correlated. All our lattice results were obtained using a non-perturbatively improved ac- 



tion [0 , and with operators renormalized on the lattice with the non-perturbative method 
of ref. [14 1, as implemented in |15, 16]. Our new result is TZ(mb)- For completeness, we also 
present some relevant 5-parameters which enter the calculation of the mixing and width 

difference 

^ 0.86(2)1™, Blf(m b ) = 0.91(3)±g : gg . (16) 



Bf s (m b ) 
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The remainder of this paper is as follows: in sec. ^| we discuss the renormalization of 
the relevant operators and the calculation of their matrix elements; the extrapolation to 
the physical points and the evaluation of the statistical and systematic errors are presented 
in sec. ^]; sec. (| contains a comparison of our results with other calculations of the same 
quantities as well as our conclusions. 



2 Calculation of the Matrix Elements 

In this section, we discuss the construction of the renormalized operators which enter the 
prediction of ATp s and describe the extraction of their matrix elements in our simulation. 
Besides the operators in eq. (^), we also need 

q s = 6*(i- 75 y^'(i- 75 y, 

Q P = b\l- l5 )s i V(l + l5 )s j , (17) 
Qp = 6 i (l-7 5 )^^(l+75)s 4 . 

The five operators in eqs. ([3D and fll7D form a complete basis necessary for the lattice 
subtractions which will be discussed later on. In ref. [pi, a new method, that allows the 



calculation of AF = 2 amplitudes without subtractions, has been proposed and feasibility 
studies are underway. If successful, it will be obviously applied also to the calculation of 

Ar Bs . 

The matrix elements which contribute to ATp s are traditionally computed in terms of 
their value in the vacuum saturation approximation (VSA), by introducing the so called B- 
parameters. The latter encode the mismatch between full QCD and VSA values. There is 
a certain freedom in defining the 5-parameters (see for example the discussion in ref. fT9||). 
For Qs and Qs, two equivalent definitions will be used in the following 

{B°\Ql(»)\B°) = \m 2 B Jl,BM, 

- 1 Uo.r+m.0.) )'^ Ssw g ■ (18) 

The first definition is the traditional one that requires, for the physical matrix elements, 
the knowledge of the quark masses; the second one may present some advantage, because 
the matrix elements are derived using physical quantities only (mg s and /b s ). The label 
(n) denotes that operators and quark masses are renormalized, in a given renormalization 
scheme (MS in our case) at the scale \i. Since the matrix element of the first operator, 
essential for B®-B° s mixing, was studied in detail in our previous paper [11], here we only 
consider the two other relevant operators, namely Qs and Qs- 
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2.1 Matrix elements from correlation functions 



Before presenting our results, let us recall the main steps necessary to extract the matrix 
elements in numerical simulations. As usual, one starts from the (Euclidean) three-point 
correlation functions: 

cf\t Pl ,t P2 ) = (o|p 5 (f,-t P2 )g,(o)p 5 (y,t Pl )|o) , (19) 

where Qi denotes any of the operators enumerated in eqs. @ and (]17D . When the Qi 
and the sources P5 (which we choose to be P5 = isj^b) are sufficiently separated in time, 
the lightest pseudoscalar-meson contribution dominates the correlation functions and the 
matrix elements can be extracted 

Cl 3 \t Pl ,t P2 ) ^ >>0 ^ge"^ (PU(a)\P) ^|e-^, (20) 

where y/Zp = |(0|P5|P)|. We take both mesons at rest and label the bare operators on the 
lattice as Qi(a), to distinguish them from their continuum counterparts. The matching to 
the continuum renormalized operators is what we discuss next. 



2.2 Operator Matching and Renormalization 

In this subsection, we describe the procedure used to get the renormalized operator Qs(n), 
relevant in the calculation of Ar# s , from the lattice bare operators. This is achieved 
through a two steps procedure: 

i) We define the subtracted operators Q' s and Q' s , obeying to the continuum Ward 
identities (up to corrections of 0(a)), as 

Q' s = Q s (a) + M<?o)<2i(a) , 

i=L,P,P 

Q' s = Qs(a)+ MoDQii*)- (21) 

i=L,P,P 

The constants A^q), which we calculate using the non-perturbative method dis- 
cussed in refs. |L5|, [Tj|, are listed in Tab. ||. Note that the mixing (^l|) is a lattice 
artifact (as a consequence of the explicit chiral symmetry breaking in the Wilson 
action) and the subtraction ensures that the resulting operators, Q' s and Q' s , have 
the same chiral properties as in the continuum. 

In principle, the mixing coefficients A^q) are functions of the bare coupling constant 
#o only. Therefore, at fixed lattice spacing, they should be independent of the scale 
at which the operators are renormalized. In practice, due to some systematic effects, 
they may depend on the renormalization scale (which corresponds to the virtuality 
of the external quark legs). This induces an uncertainty in the determination of the 
physical matrix elements which will be accounted for in the estimate of the systematic 
error. 
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Scale \x 


A L 


A P 


Ap 


A L 


Ap 


Ap 


1.9 GeV 


0.005(1) 


0.219(8) 


-0.016(8) 


-0.002(0) 


-0.094(3) 


0.007(3) 


2.7 GeV 


0.003(0) 


0.175(5) 


-0.014(2) 


-0.001(0) 


-0.075(2) 


0.005(1) 


3.8 GeV 


0.002(1) 


0.189(3) 


-0.012(2) 


-0.001(0) 


-0.081(1) 


0.003(1) 



Scale \i 


Z22 (p) 


Z23 (/*) 


Z32 (/*) 


Z33 (//) 


1.9 GeV 


0.237(13) 


-0.122(16) 


0.313(1) 


1.018(5) 


2.7 GeV 


0.282(12) 


-0.128(16) 


0.229(0) 


0.883(3) 


3.8 GeV 


0.332(12) 


-0.184(16) 


0.203(1) 


0.902(0) 



Table 2: Numerical results for the A(#q)s and i/ie matrix Z{g\,\i). They have been 
evaluated non-perturbatively in the Landau RI-MOM scheme, at (3 = = 6.2, at the 
three different scales /1 given in the table. 



CPS symmetry allows the mixing of Qs and Qs under renormalization. This is why 
we must consider both of them, although our main target is the matrix element of 
the renormalized Qs{^)- lu the second step, the operators are renormalized as: 



Qsin) \ ( ZM Z 23 (fi) \ ( Q' s 
Qs(fi) J " V Z M z 3 m J \ Q's 



(22) 



where the structure of mixing (Z 2 s 7^ Z 32 7^ 0) is the same as in the continuum. 
We compute the renormalization matrix non-perturbatively by using the method of 



ref. ID- 61, in the Landau RI-MOM renormalization scheme. The results for three values 



of the renormalization scale, // = { 1.9 GeV, 2.7 GeV, 3.8 GeV }, are given in Tab. 



2.3 Extraction of the B-parameters 

Equipped with suitably renormalized operators in the RI-MOM scheme, we proceed by 
removing the external meson propagators and sources from the correlation functions. This 
can be done in two ways. From the ratios 

cf{t Pl ,t P2 -^) (p\QM\p) = B , {) 

-lZ\Cf P {t P2 ) cf P {t Pl ) -||(0|i |P)| 2 SW ' 
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(j\QM\P) 
|l(o|i |P)l 2 



(23) 



we extract the S'-parameters (Method- J). The quality of the resulting plateaus is illustrated 
in Fig. [I]. Since the lattice renormalization constant of the axial current (Aq = Zas^q^^o) 
is ^-independent, the anomalous dimension of the parameter B' s (fi) is exactly the same as 
for(QsGu)}. 



1.5 



1.0 



Bl ([1=3.8 GeV) 



1.5 



1.0 



B' s ([1=3.8 GeV) 



******** * 



30 



35 

Time 



40 



Figure 1: The ratios defined in eq. (\29[) are shown as a function of Time = — tp 2 at 
fixed tp 1 = 16. All quantities are in lattice units. The figure refers to n q = 0.1349 and 
k q = 0.1220. 



On the other hand, if we combine the correlation functions in the ratios: 

»(3), 



C s (tp 1 ,tp 2 ; fi) 



\Z P (pi)Cpp(tp 2 ; (i) C pp (t Pl ; 



&5 (tPi > tP2 i A 1 ) 



(P\Qs(ji)\P) 

-|l(o|A(/x)|P)| 2 
(p\Qs(»)\p) _ 



(24) 

we get the standard 5-parameters (Method-IT) . The plateaus are shown in Fig. |2|. In the 
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1.0 



0.5 



B ([1=3.8 GeV) 



30 



35 

Time 



40 



Figure 2: As in Fig. |T]. for the standard B$ and B$, as specified in eq. (2^.). 



above ratios, we have used the renormalized pseudoscalar density, P*,(n) = iZp{\i)s^b. 
The value of Zp(/i) is also obtained non-perturbatively in the way described in ref. |T^| . 
Numerically (and in the chiral limit), we have Zp(n) = {0.43(2), 0.45(2), 0.50(2)} in in- 
creasing order of /i. 

The second method has the advantage that the (rather large) /i-dependence of the oper- 
ator (Qsif')) is almost cancelled by the anomalous dimension of the squared renormalized 
pseudoscalar density. On the other hand, the first method seems more convenient because 
physical amplitudes can be obtained without introducing the quark masses, which are a 
further source of theoretical uncertainty. We found, however, that B' s has a very strong 
dependence on the heavy quark mass, which prevents a reliable extrapolation to me s . Be- 
fore discussing the subtleties related to the extrapolation, we present our results for the 
heavy-light meson masses directly accessible in our simulation. 



2.4 B- parameters in the Landau RI-MOM scheme 

In this subsection we present the results for both sets of S-parameters. As in our previous 



paper |Tl[] , our study is based on a sample of 200 independent quenched gauge field con- 
figurations, generated at the coupling constant (3 = 6.2, on the volume 24 3 x 48. We use 
three values of the hopping parameter corresponding to the light quark mass (n q = 0.1344, 
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0.1349, 0.1352), and three values corresponding to the heavy quarks (kq = 0.125, 0.122, 
0.119). The first source is kept fixed at tp 1 = 16, while the second one moves along the 
temporal axis. The 4-fermion operator under study is inserted at the origin (t = 0). After 
examining the plateaus of the different ratios in eqs. (^) and fl2"4|), for every combination 
of the hopping parameters, we choose to fit in the time intervals, G [31,34] and 

tg s G [32,35]. We present results for each value of kq, with the light quark interpolated 
to the s-quark or extrapolated to the <i-quark. For a generic 5-parameter, this is obtained 
by fitting our data to the following expression: 

B(m Q ,m q ) = a® + afMp(m q ,m q ) . (25) 

This is illustrated in Fig. |], while a detailed list of results is presented in Tabs. |] and |[ 



1.3 



1.2 



1.1 



B' s (li) 



0.9 



0.05 



0.1 



0.7 



B s (\i) 



0.05 



0.1 



Figure 3: Fit in the light quark mass of Bs(fi) and B' s (fi) (at fi = 3.8 GeV) according to 
eq. $Zo]). Empty symbols denote data obtained from our simulation, whereas filled symbols 
denote quantities obtained after interpolation (extrapolation) to the strange (up-down) light 
quark mass. The heavy-quark mass corresponds to k,q 2 = 0.1220. 



At this point we have obtained the matrix elements parameterized in two ways ( B' s (n)- 
B' s (fj,) and Bs(fi)-Bs([i), respectively) in the RI-MOM scheme, at three values of the heavy 
quark masses (around the charm quark mass). The matrix elements that we need refer, 
instead, to the MS scheme and to the heavy mesons B s . Thus, to get the physical results, 
we have to discuss the scheme dependence, and the extrapolation in the heavy-meson mass. 
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n 


1.9 GeV 


2.7 GeV 


3.8 GeV 


Light quark 


q = s 


q = d 


q = s 


q = d 


q = s 


q = d 


B' s (ji- k q = 0.1250) 
5 s(a*; k Q = 0.1250) 


0.93(3) 
0.99(4) 


0.98(5) 
1.06(6) 


1.15(4) 
1.21(4) 


1.21(6) 
1.29(6) 


1.32(4) 
1.41(5) 


1.38(7) 
1.49(7) 


B' s (ji- k q = 0.1220) 
B'slfi- kq = 0.1220) 


0.83(3) 
0.87(3) 


0.85(4) 
0.91(5) 


1.02(3) 
1.06(3) 


1.05(5) 
1.10(5) 


1.17(4) 
1.23(4) 


1.20(5) 
1.28(6) 


kq = 0.1190) 
B'gQi- Kq = 0.1190) 


0.75(2) 
0.77(3) 


0.77(3) 
0.82(4) 


0.91(2) 
0.94(3) 


0.94(4) 
0.99(4) 


1.05(3) 
1.10(3) 


1.08(4) 
1.15(5) 



Table 3: Results for B s (fi) and B' s (fi) in the Landau RI-MOM scheme. The light quark 
mass is extrapolated/interpolated to d/s quarks. 



n 


1.9 GeV 


2.7 GeV 


3.8 GeV 


Light quark 


q = s 


q = d 


q = s 


q = d 


q = s 


q = d 


B S (ji- Kq = 0.1250) 
B s {fi; Kq = 0.1250) 


0.87(3) 
0.92(3) 


0.84(4) 
0.91(4) 


0.85(2) 
0.90(3) 


0.82(3) 
0.88(4) 


0.87(2) 
0.88(3) 


0.79(3) 
0.86(4) 


B s (ii; kq = 0.1220) 
Bsfa Kq = 0.1220) 


0.90(3) 
0.94(3) 


0.85(4) 
0.91(4) 


0.88(2) 
0.92(2) 


0.83(3) 
0.88(4) 


0.85(2) 
0.90(2) 


0.80(3) 
0.86(4) 


B s (n; Kq = 0.1190) 
B s (ji; Kq = 0.1190) 


0.92(2) 
0.95(3) 


0.91(3) 
0.96(4) 


0.90(2) 
0.93(2) 


0.89(3) 
0.94(3) 


0.87(2) 
0.91(2) 


0.86(3) 
0.91(3) 



Table 4: As in Tab. but for B s (fi) and B s (fi). 

3 The Physical Mixing Amplitudes 

In this section we discuss the scaling behaviour of the renormalized amplitudes (S-parameters), 
the conversion of the S-parameters from the RI-MOM to the MS scheme and the extrap- 
olation of the results in the heavy-quark mass. These points are all essential to obtain 
the final results and to estimate the systematic errors. The explicit expressions for the 
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evolution matrices and the corrections relating different schemes have been derived from 
the results of refs. p0| , pT . 



3.1 Scale dependence of the 5-parameters 

The renormalized operators obtained non-perturbatively are subject to systematic errors. 
It is thus important to check whether the renormalized matrix elements follow the scaling 
behaviour predicted by NLO perturbation theory. This is also important because we have 
finally to compute the physical amplitudes by combining our matrix elements with the 
Wilson coefficients evaluated using perturbation theory in ref. . 

The scaling behaviour of the matrix elements is governed by the following equation (we 
use the same notation as in ref. 



(QM) = W T \ii 2 ,ii x }- l {Q{ii X )) , (26) 

where the evolution operator W[//2, /^l] can be written as: 

W\pb2,nx) = M(// 2 )C/(^ 2) ^ 1 )M- 1 ( i u 1 ). (27) 

£7(/i2,/ii) is the leading order matrix and the NLO-corrections are encoded in M(/i). For 
our purpose, it is convenient to rewrite eq. (|27| ) in the following form 



W> 2 ,H =w( l x 2 )w~\fi 1 ), (28) 

where 

w(fi) = M(/i) «,(/i)-^/ 2/3 ° , (29) 

and (3q — 11 — 2n//3. The one- loop anomalous dimension is scheme independent and, in 
the basis (0), it is given by: 



7o 

The NLO contribution 



-28/3 4/3 
16/3 32/3 



m(„) = T + j£L MOM ^ (so) 

is given in terms of the matrix J^Lmomi which we write explicitly for Uj = 



^RI-MOM 



/ 170749 + 44 W 2 _ 247372 + 28 j g \ 

/ 27225 ~ 9 u & 27225 ~ 9 u & \ 

_ _6667_ , 28i 2 _ 196424 , 44 i 2 

^ 27225 ^ 9 1U 6 Z 27225 + 9 1U & Z . 



(31) 



since our lattice results are obtained in the quenched approximation. These formulae are 
sufficient for the study of the scaling behaviour of B' s {p) and B' s {p). For B s (fi) and Bs(fi), 
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since they are obtained by dividing the operator matrix elements by (P^fi)) = (0\P^(fj,)\P) , 
we also need the NLO evolution of the pseudoscalar density with the scale /i, in the RI- 
MOM scheme and with n/ = 0. This is given by 

( ^U-MOM W) = ( " 4/11 ( ^ _ m <X a {to) - *M \ ( pRI-MOM (/ii)) ( g 2) 

\a s (ni) J \ 242 7r / 

We now use the above formulae to check whether our lattice results scale as predicted 
by perturbation theory. 

In Fig. we plot the evolution of both B' s (fi) and Bs{fJ>), by normalizing them at one 
of the scales at which we have computed the renormalization matrix in eq. ( p2|) . From 
the figure, we see that B' s (1.9 GeV) falls below the other results. On the other hand, 
the evolution curves relative to B' s {ji = 2.7 GeV) and B' s {ji = 3.8 GeV) are closer to each 
other. This is to be contrasted to the situation for Bs{fi) where the scale dependence is not 
as large as for B' s (/j,), and the description of our data by the perturbative NLO anomalous 
dimension is more satisfactory, as also shown in Fig. [|. To convert our result to the MS 
scheme, as central values we choose the S-parameters obtained with the non-perturbative 
renormalization at \x = 3.8 GeV. The difference with the other results will be accounted in 
the systematic uncertainty. 



3.2 B parameters in the MS scheme at /i = m& 

In ref. [|J, the formulae in eqs. @ and (f|) were derived in the MS scheme. For this 
reason we have to convert our results from RI-MOM to MS. We have chosen to change 
renormalization scheme before extrapolating in the heavy quark mass. The change of 
scheme is obtained by using the relation 



where 



(Q MS M> 



1 + r- 



MS 



47T 



MS 



56 + 88 log 2 
- 143 + 56 log 2 



(Q 



RI-MOM 



(a0> 



-8 + 56 log2 \ 
-115 + 88 log 2 

/ 



(33) 



We note that r 
4.6 GeV using 



is independent of nj 

RI-MOM 



6|), with J™{ t!!™,, replaced by 



(Q Mb ({i = 3.8GeV)) is then evolved to fi = m b 



J-, 



MS 



V 



9561 


20723 


3025 


18150 


1811 


4997 


9075 


6050 



(34) 



/ 



m 



pole 



4.6 GeV are presented in Tab. [|. m^ ole = 4.6 GeV 



The MS .B-parameters at \i 
corresponds to the MS mass given in Tab. [1] and it is very close to the value used in 
evaluate the Wilson coefficients at NLO. 



to 
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1.5 - 



1 i.o - 



0.5 - 



2.0 3.0 4.0 

H [GeV] 



1.5 



i~ i.o 

% 



o-s L , , , , , , , , , , i Z 

2.0 3.0 4.0 

H [GeV] 



Figure 4: NLO evolution of the parameters B' s (fi) and Bg{fi) in the RI-MOM scheme. 
The curves are obtained by starting the evolution from the lattice results (filled circles) at 
\i = 1.9, 2.7 and 3.8 GeV, respectively. 
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0.1250 


0.1220 


0.1190 


Light quark 


q = s 


q = d 


q = s 


q = d 


q = s 


q = d 


mp [GeV] 


1.85(7) 


1.75(8) 


2.11(9) 


2.02(9) 


2.38(10) 


2.26(11) 


B' s (m P ,m b ) 
B f s (m P ,m b ) 


1.52(5) 
2.18(8) 


1.60(8) 
2.31(12) 


1.35(4) 
1.93(6) 


1.39(7) 
1.99(10) 


1.21(3) 
1.72(5) 


1.25(5) 
1.79(8) 


B s (m P ,m b ) 
B s (m P ,m b ) 


0.73(2) 
1.05(3) 


0.71(3) 
1.03(5) 


0.76(2) 
1.08(3) 


0.72(3) 
1.03(5) 


0.78(2) 
1.10(3) 


0.76(3) 
1.09(5) 



Table 5: B '-parameters in the MS scheme at /x = 4.6 GeV. 
3.3 Extrapolation to the B s meson 

From Tabs. [3], |3] and we observe that the dependence of B' s on both the renormalization 
scale and the heavy-quark (meson) mass is much more pronounced than in the case of the 
parameter B s (see also Figs. § and |). In particular, because of the large mass dependence, 
it is more difficult to extrapolate B' s to the physical point. This is related to the fact that 
in eqs. (|T8|) the ratio (m^/mi,) 2 , distinguishing Bs from B' s , still varies very rapidly in the 
mass range (around the charm) covered in our simulation. A similar problem is found if 
one tries to extract m b from the heavy-meson spectrum, computed with fully propagating 
quarks, by extrapolating in the heavy-quark mass. For this reason, so far, m b has been 
computed on the lattice only by using the Heavy Quark Effective Theory (HQET) |22| 
or Non-Relativistic QCD (NRQCD) ||23|| . Thus, although in principle we would prefer 
B' s because it allows the evaluation of the matrix element without using the quark mass, 
in practice our best results are those obtained from Bs- For the same reasons, we were 
unable to extract 1Z{m b ) = (Qs)/ (Ql) directly from the ratio of the matrix elements. The 
strongest dependence of B' s on the quark mass finds its explanation in the framework of 
the HQET, by studying the leading and subleading contributions in the l/m b expansion 
to the B parameters Our estimates of (Ar^/r^J are then obtained with 7Z(m b ) 

computed by using Bs, Bb s and m b from Tab. [1], as shown below. For completeness we 
also present the results for B' s . Hopefully, when smaller values of the lattice spacing, and 
hence larger values of the heavy quark mass, will be accessible, an accurate determination 
of the physical matrix element will be provided by B' s . 

In order to make the extrapolation in the heavy-quark mass from the region where we 
have data to m b , we rely on the scaling laws of the HQET. The S-parameters in Tab. |^ 
have been obtained using the anomalous dimension matrix of the massless theory. This is 
the appropriate procedure since, on our data, /i is larger than the heavy quark mass, m Q , 
used in our simulation. In order to use the HQET scaling laws, we have to evolve first to a 
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scale smaller than the quark mass, at tuq fixed, and then to study the scaling at fixed fi as 
a function of the heavy quark mass. In practice, this is achieved, at LO in the anomalous 
dimensions, by introducing the following quantities: 

$>(m Ps ,m b ) = s ) Ps) B(m Ps ,m b ) 
\a s {m Bs )J 

$'(mp 3 ,m 6 ) = ( °B'(m Pa ,m b ) (35) 



a s (m B 

where we have introduced the matrices 7 = 70 — 7o + 2(tp — jp), and 7' = 70 — 7o + 27a- 
B = (^Bs, BsJ and similarly for B' . 70 is the HQET anomalous dimension matrix which 



at leading order is given by 



7, ■ -- \ \ \ ) ■ (36) 



whereas 7p = 7a = — 4 and the LO anomalous dimension of the pseudoscalar density in 
full QCD is given by 7p = —8. 

The quantity Q(m Pa ,m b ) scales with the heavy quark (heavy meson) as: 

$(m Ps ,m fe ) = a + — + ... (37) 
m Ps 

The extrapolations of B^lrrib) and B'^lrrib) to m Pa are shown in Fig. |^. The results 
are: 

Bf*(m h ) = 0.86(2) , 5f^(m 6 ) = 1.25(3) 

B' s m (m b ) = 0.63(3) , B' s m (m b ) = 0.91(5) . (38) 
As for the systematic uncertainties, we account for the following two: 

• if, instead of converting to the MS scheme before extrapolating to m Ps , we extrapolate 
our RI-MOM results (obtained at /1 = 3.8 GeV) in the heavy quark by using the same 
scaling laws of eqs. (p5|) , then evolve to /i = m b , and convert to the MS scheme, the 
results are 0: 

£fS(m 6 ) = 0.88(2) , 5^(m fc ) = 1.27(3) 

B' s m (m b ) = 0.62(3) , B' s m (m b ) = 0.90(5) . (39) 

• if we start from a fi lower than 3.8 GeV, for example [i = 1.9 GeV, then Bs(m b ) and 
Bs(m b ) change by —3%, while B' s {m b ) and B' s {m b ) drop by —6%. 



3 If the last conversion from the RI — MOM to the MS scheme is made by using rif = 4 the results 
remain remarkably stable (they decrease by about 1% ). 
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0.75 
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♦ 3>'(m p ,m b ) 
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1.00 1.25 1.50 

l/m„ 



Figure 5: Extrapolation of §(mp,m b ) and $'(mp,m(,) in 1/rrip to the B s -meson mass, mp 
is given in lattice units. 



After combining these uncertainties, we arrive at our final results: 

B^(m b ) = 0.86(2)^;^ , = 1.25(3) +0 - 02 



0.05 > 



B' s MS (m b ) = 0.63(3) +0 - 00 



0.04 > 



B' s MS (m b ) = 0.91(5)18:" 



(40) 
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From these numbers, we observe that there is a discrepancy between the value of 
Bf(m & ) and B' s m (m b ). Using B' s m (m b ) ~ {m B Jm b fB s ls {m b ), w i t h (m B Jm b ) 2 = 1.6 
from Tab. |I|, we get B' s (m b ) ~ 1.37, which is about twice the value of B' s MS (rrib) in (JTO|). 

In order to obtain the physical results we have used B^ s (rrib), for which the extrapola- 
tion to the B s meson is much smoother. In this case, we are obliged to use the quark mass 
to derive the matrix element needed to obtain TZ(mb), which is the relevant quantity for 
(ATbJTbs), as explained in the introduction. This is not a major concern, however, since 
nib is known with a very tiny error (see Tab. |l|). 

The best way to minimize the statistical uncertainty is to compute on the same set of 
configurations 

Sr^ = °- 95 ( 3 )-o:o2 , (4i) 

which, when combined with the numbers from the Tab. [l|, gives the wanted quantity 

5/ m Bs V B s {m b ) n rtO/oX+a00 



Jl(m b ) = -- — -, — x Ds , — r tt^V-^t = -0.93(3) + ^^ . (42) 
1 b> 8 \m b (m b ) + m s (m b ) J B Bs (m b ) 1 ; -°' 01 1 ; 

4 Conclusion 

In this section, we compare our results with previous lattice studies and draw our conclu- 
sion. 



The first lattice calculation of Bs(fJ-) was performed in ref. [2"5fl, using the (unimproved) 
Wilson fermion action. Their 5-parameters were obtained in a modified MS NDR scheme, 
at the scale /i = 2.33 GeV, at a heavy-quark mass corresponding approximately to the D s 
mass. By translating their result to the MS of ref. §, one finds E!p(2.33 GeV) = 0.81(1) 
and B S AS (2.33 GeV) = 0.86(1). To compare with these numbers, we have used our results 
at k,q = 0.1250, interpolated to the strange light-quark mass and then evolved down to 
fi = 2.33 GeV. We get 5fS(2.33 GeV) = 0.80(2) and 5^(2.33 GeV)^ 1.05(2). Bf* 
is hence in very good agreement with the value of [^|, whereas our B^ is 20% larger. 
We note, however, that the values 5^(2.33 GeV) = 0.81 and 5^(2.33 GeV)_= 0.86 
correspond (only evolution but no extrapolation in the heavy quark mass) to Bf s (m b = 
0.75 which is 14% smaller than our number. 

In refs. [17 , |26[| , NRQCD has been used to compute B' s (m b ). A comparison with them 



is interesting, because effective theories have different systematic errors with respect to 
the approach followed in the present study. Their latest result (which we convert to Bs) 
reads i?^ s (mb) = 0.78(2)(10), which is in fair agreement with ours. Moreover their ratio 
Bg(m b )/Bg(m b ) ~ 0.78/0.85 = 0.92 is very close to ours. This implies that, in spite of 
the difference in the single 5-parameters, the same physical answer for (ATb s /Tb s ) is ob- 
tained from eq. (§). In |L7], |26|, in order to get the width difference, they adopted, however, 
another expression (eq. (4) of [Pq| ), which uses different inputs, namely the experimental 
inclusive semileptonic branching ratio and the theoretical determination of the decay con- 
stant f Bs = 245±30 MeV. Their result 0, (Ar B ,/r B ,) = (10.7±2.6± 1.4± 1.7) x 10~ 2 is 
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about a factor of two larger than ours. One may wonder whether the difference is due to a 
different evaluation of the 1/m corrections, which are so important in this game. From our 
calculation, we find that this contribution to (ATbJTbJ is about —8.6 x 10~ 2 , identical 
to the value used in |L7] and (eq. (9) of f26f ). Besides the fact that they use a for- 



mula which involves the inclusive semileptonic branching ratio, the main difference stems, 
instead, from the use of a very large value of /b s (and to a lesser extent from the use of 
m b = 4.8 GeV instead of 4.6 GeV). We do not find it justified to use the unquenched value 
of f'B s , combined with B parameters computed in the quenched approximation. What 
really matters is the combination of these quantities in the matrix elements, and we do 
not know how much the B parameters change in the unquenched case. This is the reason 
why we prefer eq. (§), which does not require fB s , but only lZ(rrib) (which is essentially the 
same for us and in |T7[]) and £, which is known to remain almost the same in the quenched 
and unquenched case |10[ . 



On the lattice, physical quantities relevant in heavy-quark physics can be computed 
following two main routes, either by extrapolating the results in the heavy-quark mass from 
a region around the charm mass or by using some effective theory (HQET or NRQCD). 
The two approaches have different systematics and in many cases lead to results which 
are barely compatible. The B s width difference is particularly lucky, in this respect, since 



TMjrib) in our calculation and in ref. |T7| are in excellent agreement and lead to the same 
value of (ATb s /Tb s ) if eq. (|8|) is used. It is not surprising that ref. |HJ predicts a much 
larger value for (ATbJTbJ, since they use a very large value of /# s , which is not needed 
in eq. @. We find that, in order to reduce the present error, a better determination of 
the 1/m correction, although rather hard, is very important. Obviously, calculations with 
larger heavy-quark masses and without quenching are also demanded for a better control 
of the remaining systematic errors. 
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